Paired state in an integrable spin-1 boson model 
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An exactly solvable model describing the low density limit of the spin-1 bosons in a one- 
dimensional optical lattice is proposed. The exact Bethe ansatz solution shows that the low energy 
physics of this system is described by a quantum liquid of spin singlet bound pairs. Motivated by 
the exact results, a mean-field approach to the corresponding three-dimensional system is carried 
out. Condensation of singlet pairs and coexistence with ordinary Bose-Einstein condensation are 
predicted. 
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INTRODUCTION 

Study on the trapped cold atoms opens the door for 
finding new matter states which are usually unknown or 
even do not exist in nature. Experimentally, the cold 
atom gas can be realized by means of either magnetic or 
optical traps. With Feshbach resonance, the scattering 
length and thus the couplings among atoms can be ma- 
nipulated experimentally. In addition, with laser beams, 
one can confine particles in valleys of periodic potential 
of the optical lattice. These experimental tools provide 
a platform to study quite clean and controllable "artifi- 
cial condensed matter systems" . Moreover, particles with 
higher inner degrees of freedom (hyperfine spin), which 
usually do not exist in conventional condensed matters, 
can be prepared by catching several hyperfine sublevels of 
atoms. Compared to spinless Bose gases, the low-energy 
physics of these systems such as the spin dynamics [l|,|2|,|3( 
is much richer and may show fascinating macroscopic 
quantum phenomena. For example, the multi-component 
Bose-Einstein condensation (BEC) is realized in 87 Kb [4| 
and 23 Na @, @, 0] gases with optical traps. Both 87 Rb 
and 23 TV a atoms have a hyperfine spin F — 1. The inter- 
action among 87 Rb atoms is ferromagnetic, which leads 
to a spin-polarized (ferromagnetic) ground state, while 
the spin exchange interaction among the 23 N a atoms is 
antifcrromagnetic, leaving the ground state a spin sin- 
glet condensate. In an optical lattice, the Mott phase of 
F = 1 cold atoms may exhibit rich magnetic structures. 
Nematic singlet [8[ or dimerized [9( ground state has been 
proposed. Nevertheless, study on spinor cold atoms is 
still young and a quite interesting issue 1(J U, 12, HH 
in modern many body physics. 

In this Letter, we propose an exactly solvable model 
for F = 1 bosonic cold atoms. The Bethe ansatz so- 
lution exactly shows that atoms may form spin singlet 
pairs with a finite energy gap and the low-energy physics 
is described by a quantum liquid of spin singlet atom 
pairs. Based on the exact solution for the ID model, an 
appropriate mean-field theory is proposed to study the 
corresponding 3D systems. BCS-like pair condensation 



and coexistence with ordinary BEC are found in the 3D 
model. 



THE MODEL 

In an optical lattice, it has been pro posed that the 
following boson Hubbard model 1J, [l5| well describes 
the low-energy physics of the spinor bosons: 

<ij>,s i 



H 



where a is (at s ) is the creation (annihilation) operator 
of atoms on site i with spin index s, rij and are the 
particle number and spin operators, respectively; \i is the 
chemical potential. 

Recently, tremendous experimental and theoretical 
progress has been achieved in realization of one- 
dimensional (ID) cold atom systems [3, 17, [HI, 
|2(I l2l| and 2D systems [22] . The Mott phase diagram 
of F = 1 bosons in an optical lattice has been given 
in ref . fl3j . In the metallic phase, it is known that the 
ID spinless bosonic atom_gases are well described by 
Lieb-Liniger model [2S" 



24| and several physical prop- 



erties based on Lieb-Liniger 's exact results have been 
derived [2 5l 1 2 61]. However, results on ID cold atoms with 
internal degrees of freedom in the metallic phase are still 
rare. Generally speaking, a ID exactly solvable model 
not only gives the best understanding for the correspond- 
ing universal class, but also provides some useful clues for 
understanding three-dimensional (3D) systems. 

In this paper, instead of studying model (1), we con- 
sider the following ID Hamiltonian: 



N 



ff= "E^2 + Et c ° + C2Si ' S il S ( x i ~ x i)> ( 2 ) 

i— 1 1 i<j 

where Si is the spin operator with ^-components s =f 



,0,|; c 



(.9o + 232)/3, c 2 = {g 2 - 5o)/3, gs 
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4irh 2 ls/Mi,, Mb is the mass of boson and Is is the s-wave 
scattering length in the total spin 5 channel (HQ. I n the 
second quantization form, we define the particle creation 
(annihilation) operators as a\(x) (a s (x)). Obviously, the 
model (2) is just the low density limit of the boson Hub- 
bard model (1). We note that two-particle scattering 
processes keep the conservation of the total spin 5 and 
therefore the model possesses an SU(2) invariance. Non- 
trivial scattering occurs only in the 5 = and S = 2 
channels. In the 5 = 1 channel the wave- function is an- 
tisymmetric by exchanging two particles and the delta- 
function interaction is irrelevant. Especially in the 5 = 
channel, a special scattering process 

a\(x) +a\(x) -> 2a\(x) (3) 

occurs, which makes the total particle number of an in- 
dividual spin component is no longer a good quantum 
number and breaks the SU(3) invariance. It is easy to 
verify that the present model has the following conserved 
quantities: 

N - e/^m*. 

S 

S z = f[a\(x)a^x)-a\(x)a L (x)]dx, (4) 

where N and S z are the total particle number opera- 
tor and z-component of the total spin operator, respec- 
tively. Because of the SU(2) invariance of the Hamilto- 
nian, there are also two other good quantum numbers: 

5 + = V2 J [a^(x)ao(x) + aQ(x)ai(x)]dx, 

5~ = V2 J [aj ) (x)a^(x) + a^(x)ao(x)]dx. (5) 

S z and form the generators of the SU(2) algebra. 
These three spin operators, combined with the five spin 
quadrupole operators 



Qq = j {a^(x)a^(x) + a^(x)ai(x) ~ 2ali(x)ao(x)]dx, 
Q 2 = J[a^(x)ai(x) + a^(x)a^(x)]dx, 
Q xy = -i / [a^(x)ai(x) - a^(x)a^(x)]dx, (6) 



Qxz = A= J[a\(x)a Q (x) - al(x)ai(x) +h.c.]dx, 



Q yz = J [a\(x)a (x) - al(x)ai(x) - h.c.]dx, 

form the basic representation of the SU(3) algebra. 



BETHE ANSATZ SOLUTION 

The pioneer work on the integrable models with inter- 
nal degrees of freedom was done by Yang 27 , and fol- 



lowed by Sutherland [29]. There are two integrable lines 
for the model (2). The first is the C2 = case, i.e., SU(3)- 
invariant case, which has been solved by Sutherland [29(. 
The second integrable line is cq = C2, which has never 
been studied before and is the main target of the present 
work. In the framework of coordinate Bethe ansatz, the 
wave function of the system described by a set of quasi- 
momenta {kj} can be written as[27l [28j 

*(xisi, • • • ,x N s N ) = ^2 ®( x Qi < ■ ■ < xq n ) 
Q,P 

xA Sl ... SN (Q,Py^ k ^, (7) 

where Q = (Qi, ■ ■ ■ ,Qn) and P = (Pi, • • • , Pjy) are the 
permutations of the integers 1, • • • ,N, 9(xq 1 < ■ ■ ■ < 

X Q N ) = 0( x Qn -xQn-i) ■ ■■0( X Q2 ~ x Qi) and G{x-y) is 
the step function. The wave function is symmetric under 
permutating both the coordinates and the spins of two 
atoms. The wave function is continuous but its derivative 
jumps when two atoms touch. With the standard coor- 
dinate Bethe ansatz procedure, we obtain the two-body 
scattering matrix for cq = C2 = c as 



ic 



2ic p2 

2ic ij ' 



(8) 



where S = 0,1,2 is the spin projection operator 

onto the state of total spin 5. The scattering matrix 
satisfies the Yang-Baxter equation [27L Eil. liij] 

5i 2 (fci - fc 2 )5i 3 (fci - fc 3 ) 5 23 (fc2 - h) 
= 5 23 (fc 2 - fc 3 )Si3(*i ~ k 3 )S 12 (k! - k 2 ), (9) 

which ensures the integrability of the model ^ at cq = 
c 2 = c. With the periodic boundary conditions of the 
wave function, we obtain the following eigenvalue equa- 
tions 

SjNSjN-i ■ ■ ■ Sjj+iSjj-i ■ ■ ■ 5jie lfc:,i £;o = Co, (10) 

where £0 is the amplitude of initial state wave function. 
We follow the algebraic Bethe ansatz method developed 
in [30L [3lL [3^ | to solve the above eigenvalue problem. In 
fact, the 5-matrix of the present model has the same 
structure to that of the i?-operator of the Takhtajan- 
Babujian model[3C' 



31j . In such a sense, the spin dy- 



namics of our model keeps some similarity to that of the 
Takhtajan-Babujian spin chain. Firstly, we define the 
monodromy matrix as 

T,(X) = SojSoNSoN-i ■■■ Soj+iSoj-i ■■■ 5oi 
Ai(\) Si (A) B 2 (X)\ 
Ci(A) A 2 (X) B 3 (X) , (11) 
C 2 (A) C 3 (A) A 3 (A) J 

where 5o; = 5o/(A — ki). The eigenvalue problem (10) is 
therefore reduced to 



tr T l (k j )e lk > L Zo = Z . 



(12) 
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The monodromy matrix satisfies the Yang-Baxter rela- 
tion 



on the vacuum state gives nonzero values and thus can 
be regarded as generating operator of eigenstates 



S 12 (A - u)Ti(X)T 2 (u) = T 2 (u)T 1 (X)S 12 (X - u). (13) 
Further, we define an auxiliary monodromy matrix as 

rp/\\ nas nas nas nas nas 

1 \ A ) — D 0j D 0N ' ' ' "-'Oj+l'-'Oj-l'-'Ol 



with 



SS'(A) 



( A{\) B{\) 
\C{\) 25(A) 



A — ki — ihc— ic&o • S; 
A — ki + j~c 



(14) 



(15) 



The monodromy matrices (fTT|) and (fT4|) satisfy the Yang- 
Baxter relations 

S%(\ - u)Ti(A)T 2 (u) = T 2 (u)Ti(A)S£(A - u), 
S^(X - M )Ti(A)T 2 ( M ) = r 2 ( M )Ti(A)S£f (A - m), (16) 

with S??(X) = (A - ic) -1 (A - ic/2 - ica x ■ a 2 /2). From 
Eq. (|16j) we obtain the following commutation relations 

A — u + ilc 

Ai(A)2?(u) = 



A — u — z±c 



A — u — i^c 



5-Bi(A)i4(«), 



(17) 



A 2 (X)B(u) = 



(X-u + i|c)(A - it - if c) 



(A — u + isc)(A — u — ihc) 



-j-^-B(u)A 2 (X) 



lV ~ 2c . x Bi(X)D(u) iV ^,, B 3 (X)A(u) 
X — u — i^c X — u + i^c 

2ic 

(A — u + i^c)(A — u — i^c) 2 ^ ^ ' 



(18) 



A — ?; — i-c 
A 3 (X)B(u) = - — \ B(u)A 3 (X) 



X — u + i-ljC 



iy/2c 



X — u + i^c 



—B 3 (X)D(u). 



(19) 



Meanwhile, the commutation relations of A(A), D(X) and 
B(X) read 

A — ii -4- i c i c 
A{X)B(u) = - B(tt)A(A) - t B(A)A(u)(20) 



A- 



A- 



A — 7/ — - ? f i c 

D{X)B(u) = — B(u)D(X) + - B(X)D(u%21) 



X — u 



X — u 



The vacuum state of the system is defined as |f2) = 
| |)i ® ' ' ' ® I T)jv- It is a common eigenstate of Ai(A), 
A 2 (A), A3 (A), A(A) and D(X). The element C(A) acting 
on the vacuum state gives zero. The element 23(A) acting 



\y) = B( Ul )---B(u M )\n). 



(22) 



tr T(kj) = J2n=iA n (kj) acting on the assumed Bcthc 
states (|22|) gives two kinds of terms, i.e., wanted and 
unwanted terms. Putting the unwanted terms to be zero 
we readily obtain the following Bethe ansatz equations 
for the rapidities {kj}, 



n kj - ki + 2ic i r kj - Ag-ic ) . ) 



N 

n 



, kj k\ 2ic A Q + ic 

l=l,l^j J a=l J 

M 



A a - h - ic 



n 



A Q — Ap — ic 



\* A a -ki+ ic A Q - A/3 + ic 



(24) 



where j, Z = 1, • ■ ■ , TV, a, (3 — 1, • • • , M, M is the number 
of flipped spins and {A Q } is the set of the spin rapidi- 
ties. The corresponding eigenvalue of the Hamiltonian 
(2) reads 



JV 



(25) 



THERMODYNAMIC LIMIT 

Above we have confined the particles in a finite ID 
box with length L. Based on the solutions of the Bethe 
ansatz equations, we can study the ground state and low- 
temperature properties of the system in the thermody- 
namic limit L — ► 00, N/L — ► n. The solutions of the 
Bethe ansatz equations are a little bit complicated. Be- 
sides real solutions of {kj}, {A Q }, Eqs. (|23ti24|) have also 
complex solutions for both c > and c < 0, which are 
usually called as string solutions. For c > 0, attrac- 
tive interaction only occurs in the 5 = channel. That 
means particles may form spin singlet bound pairs. Gen- 
erally, the complex solutions are determined by the poles 
or zeros of the Bethe ansatz equations in the thermo- 
dynamic limit. For example, if some kj in the upper 
complex plane, the left side of Eq. (j2"3"|) tends to zero 
when L — > 00. Correspondingly, there must exist a A a 
satisfying kj — A a — ic — > 0. Furthermore, from Eq. ([24]) 
we learn that there is another A^ with A a — A^ + ic — > 0. 
For the complex conjugate invariance of the equations, 
we obtain the simple conjugate kj-pair solutions 



kj = Kj +ic/2 + o{e- SL ), 
kj = Kj - ic/2 + o(e- SL ^ 



associated with A 2-strings 



Aj = Kj + ic/2 + o(e~ s ' L ), 



Kj - ic/2 + o(e- s ' L ), 



(26) 



(27) 
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where Kj is a real parameter, 6 and 5' are some positive 
constants. We studied the 3 and 4-particle cases and 
verified that Eq. (j2U|) describes the only possible bound 
state in the charge sector. In fact, no more than two 
atoms can form a bound state because of the symmetry 
constraint of the wave functions. In the thermodynamic 
limit, each bound pair contributes bound energy A = 
c 2 /2. Therefore, the low energy physics of the present 
system must be described by a quantum liquid of these 
bound pairs. In the ground state, all particles form such 
kind of bound pairs (Even N is supposed. For odd N, 
there is a single unpaired particle and the ground state 
is 3- fold degenerate). Substituting these 2-string ansatz 
into Eqs. (|23H24j) and taking logarithm, we arrive at one 
set of reduced Bethe ansatz equations 



JV 



K j L = *I j - ]T 



arctan 



arctan 



Ki - K, 



arctan 



2(Kj - K t 
3c 

2(Kj - K t 



, (28) 



where Ij is integer (half integer) for N/2 odd (even). The 
ground state corresponds to a sequence of consecutive 
7/s around zero symmetrically. In the thermodynamic 
limit, define po(Kj) = L^ 1 (Ij+i — Ij) / '{Kj+\ — Kj)) as the 
density of flipped spins. Taking derivative of Eq. ([28]) . we 
obtain that the density distribution po {K) in the ground 
state satisfies the following integral equation 



1 1 f c 
p (K) = - + - \ 

7T 7T J_, 



6|c| 



9c 2 + A{K - K') 2 
21c! 



Po (K')dK',{29) 



c 2 + (K - K') 2 Ac 2 + {K - K') 2 
where the pseudo Fermi point Q is determined by 



n = 2 / p (K)dK. 

J-Q 

The density of the ground state energy reads 



Q 



2K 2 -- )p (K)dK. 



(30) 



(31) 



Obviously, the ground state is a global spin singlet with 
S z = N - M = 0. However, it is not an 577(3) sin- 
glet state. In the insulator phase of the boson-Hubbard 
model, a s pin nematic state or a spin quadrupole polar- 
ized state |33T | has been obtained. In our case, the spin 
part of the wave function of each bound pair takes the 
form (| t)4) + l X ^ T) — |0, 0) )/ a/3- It can be easily deduced 
that the expectation values of the quadrupole momenta 
per unit length are (Qo) = —n/2 and (Q a ) — for a^O. 
There is a finite energy gap A = c 2 /2 for the spin excita- 
tions. The only basic gapless excitation is in the charge 
sector. This can be realized by either digging a hole in 



the K pseudo Fermi sea or putting a particle above the 
pseudo Fermi point (Note we treat an atom pair as a sin- 
gle particle here). The excitation energy e(K) of a hole 
or a particle with quasi momentum P{K) = ■nI(K)/L 
satisfies 



e(K) = 2{K 2 



') + 



1 



Q r 

-Q 
2\c 



6 c 



9c 2 + 4(K - K') 2 

e(K')dK'. (32) 



c 2 + (K - K') 2 Ac 2 + (K- K') 2 

Other gapless excitations such as the particle-hole and 
current excitations can be expressed as the superposition 
of a single particle and single hole excitations. The Fermi 
velocity is 



e '(Q) 
npo(Q)' 



(33) 



At low temperatures T -C A, the spin degrees of freedom 
are frozen completely Thus the low temperature physics 
is almost the same to that of the Lieb-Liniger model. As 
a Luttingcr liquid, its low-temperature specific heat and 
susceptibility behave as 



C(T) = — T + o{T 2 ), X (T) 
3vp 



(34) 



where we have taken the Boltzmann constant ks as our 
unit. 

The general excited states are characterized by a set of 
real {kj}; a set of k — A pairs described by Eqs. (|26H27|) 
and A n-strings taking the form of A^ = +i(n + 1 — 
o(e~ SL ), where j — 1, • 



2j>/2- 
with n 



, n and a = 1 , • • • , M n 
= 1,2, ••■ and M = ^ n nM n . Denote a', p 
and a n as the densities of bound pairs, real rapidities 
and A n-strings in the thermodynamic limit, respectively, 
and o' h , p h and <j\ as the correspon ding hole densities. 



By minimizing the Gibbs free energy |34|, we obtain the 
following coupled nonlinear integral equations 

r 2 

hxrf = 2T-\k 2 - — - (j,) - (a s - oi) * ln(l 



— (a 6 + a 4 - a 2 ) * ln(l + r/ 1 ), 
InC = T^ik 2 -(j,-2h)- (a s - oi) * ln(l 

+ lnr?i - (04 + a 2 + a a ) * ln(l + ^ _1 ), 
hxrn = G* [ln(l + 7 ?2 ) + ln(l+r 1 )], 



(35) 



In r} n = G* 

,. hi7j„ 
lim 



[ln(l 

_ h 

- 5;, 



ln(l + 7?„_i)],n = 2,3, 



+ {Axf)H = a' h /a',i = 
c 1 sech(2irx / c) , f * g — 



where a n (x) = 4n|c|/[7r((nc) 2 

P h /P, Vn = <Jn/Vn, G{x) = 

J f(x—y)g(y)dy, and h is the external magnetic field. For 
T = and h = 0, it is easily to deduce that p = a n = 
and a' = po . This also confirms that the previously given 
ground state is the correct ground state. 
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For c < 0, the interaction in S = 2 channel is attractive 
while it is repulsive in the S = channel. From the 
Bethe ansatz equations we learn that the ground state is 
a incompressible ferromagnetic state described by an N 
string 



4.0 



k 3 = ic(N+ 1 - 2j), j = l,2---,iV. 



CORRESPONDING 3D MODEL 



(36) 



Now let us turn to the 3D case. An obvious fact is 
that two kinds of condensation may occur in the corre- 
sponding 3D systems with attractive interaction in the 
5 = channel. One is the conventional BEC and the 
other is the BCS like pair condensation as indicated by 
the ID exact result. An interesting question arises: Is 
there any BCS-BEC crossover or BCS-BEC coexistence? 
To answer this question, we consider the following Hamil- 
tonian 

H = -J2 J a|(r)V 2 a s (r)dr - v J p^(r)p(r)dr, (37) 

where p*(r) = [a|(r)aj(r) + aj(r)a|(r) - aJ(r)aJ(r)]/\/3 
and v is a positive coupling constant. For simplicity, 
repulsive interaction in the S = 2 channel is omitted since 
it is irrelevant to the pair condensation. Motivated from 
the ID exact result, we introduce the order parameter of 
pair condensation as 



O 



p(r)dr 



(38) 



where V and 



denote the volume and the thermo- 



dynamic average, respectively. By using BCS-like mean- 
field approximation, we linearize (|37p as 



H 



E E £ ( k H( k K( k )-^ 

-aj(k)aj(-k) + h.c. } + VvO 



E 



4(k)4(-k) 



(39) 



where e(k) = k 2 , a =T,I and a means the spin flipped 
state. With the following Bogoliubov transformations 

6+(k)=«(k)4(k)+«(k)a 9 (-k) > 

6j(k) = u'(k)4(k) + i/(k)ao(-k), (40) 

where u 2 = u' 2 = (g + l)/2, v 2 = v' 2 = (g - l)/2, 
uV = -uv = g/f, g = fljf^l and / = V3k 2 /(vO), 
the Hamiltonian (f3"9"| can be diagonalized. The order 
parameter O and the chemical potential fi are determined 
by the following self-consistent equations 



1 

1 

4^2 



coth f^-de, 
E 2 ' 



£F We 
2 



_^ coth ^_ , ) (/( . N |, 




v=4 


PC 


v=4 


BEC 


- - v=2 


PC 


v=2 


BEC 


• v=0 


BEC 



0.6 



FIG. 1: The transiton temperatures T£ and T' versus density 
of particles n for different interacting strength v. The solid 
and dash-dot lines indicate T^, while the dashed and dotted 
lines indicate T^. The asterisk line is for the ideal gas 
(v = 0). Tl > T h c for any positive v. 



where ep is thc 
an optical lattice, 



energy 
E = 
•)dr) 



cutoff 



or band width in 



H) 2 -4v 2 2 /3, n = 
[r)dr) T is the density of particles, (5 = 
1/T. The critical temperature Tj? for pair condensation 
is determined by Eq. |4T|) with 0\t^t? — 0. Inter- 
estingly, besides the condensation of atom pairs, ordi- 
nary BEC also occurs at low temperature T <T\ when 
fi = 2vO/V3. Thc numerical solutions of T h c and T% for 
v = 2,4 are depicted in Fig. [TJ It is shown that both 
and TJ? increase with increasing interacting strength 
v and the density of particles n. Meanwhile, for a fixed 
v, TP is always larger than T c fc . When v -> 0, -> T b c . 
Below T c & , coexistence of pair condensation and BEC oc- 
curs. However, there is no BCS-BEC crossover which 
usually occurs in fermion gases. 



CONCLUSION 

In conclusion, we propose an exactly solvable model 
describing the low density limit of the spin-1 bosons in a 
ID optical lattice. Based on the exact result, a mean-field 
approach for the corresponding 3D model is introduced. 
A new matter phase, i.e., the pair condensate and coex- 
istence with ordinary BEC are predicted. We expect this 
new matter state could be realized in experiments. 

This work is supported by NSFC under grant 
No.10474125, No.10574150 and the 973-project under 
grant No.2006CB921300. 
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